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A simple extension of existing models for protein crystallisation is described, in which salt ions 
and charge neutrality are explicitly incorporated. This provides a straightforward explanation for 
the shape of protein crystallisation boundaries, the associated scaling properties seen for lysozyme, 
and can also explain much of the salt dependence of the second virial coefficient. The analysis has 
wider implications for the use of pair potentials to understand protein crystallisation. 



I. INTRODUCTION 



Protein crystallisation is of great practical importance, 
both as a purification method and since high quality crys- 
tals are needed for X-ray diffraction work ||^, ^. Under 
typical conditions where crystals are obtained, protein- 
protein interactions appear to be characterised by hard 
core repulsions with short-ranged attractions H, ||. For 
instance, Rosenbaum et al |^ successfully collapsed the 
crystallisation boundaries for a number of proteins onto 
the adhesive hard sphere (AHS) model 0, || by match- 
ing the second virial coefhcient B2- Whilst the AHS 
model is frequently used to understand protein crystalli- 
sation, it is also possible that crystallisation is an ener- 
getic ordering transition driven by highly directional in- 
teractions. A model of hard spheres with sticky patches 
was introduced by Sear to capture this important possi- 
bility §. 

In this paper I examine the properties of a simple 
model for protein crystallisation which takes into account 
the effects of charge in a very elementary way. Such a 
model suggests that some generic aspects of protein crys- 
tallisation can be explained as a straightforward conse- 
quence of charge neutrality, at least for lysozyme which 
readily undergoes crystallisation. These various aspects 
are: (i) the overall shape of the crystallisation bound- 
ary, (ii) a scaling collapse of the crystallisation boundary 
noticed by Poon et al |l0| when plotted as a function 
of Cs/Q'^, where Cg is the NaCl concentration and Q the 
protein charge; (iii) a similar scaling of the lysozyme sol- 
ubility data in Guo et al [ pT] | ; and (iv) a similar scaling 
for B2 noticed by Egelhaaf and Poon ||l2[ . 

The appearance of Cs/Q^ in all these is very remini- 
cent of the Donnan membrane equilibrium |T3], where 
one finds a contribution to B2 equal to jA^Cs provided 
one is in the high-salt limit Cs S> Qcp where Cp is the pro- 
tein concentration. The effect is well known in polyelec- 
trolyte theory . This result is usually derived by 

supposing ideal solution behaviour and imposing electri- 
cal neutrality on the two sides of a semi-permeable mem- 
brane. In the present case, I argue there is an analagous 
Donnan equilibrium between protein crystals and a dilute 
protein solution, with the crystal-solution interface play- 
ing the role of the membrane. The appearance of Cg/Q"^ 
as a scaling variable reffects the fact that the solution (not 



necessarily the crystal) is in the high-salt limit, as will be 
explained below. Donnan's theory appears to have little 
to do with the more conventional McMillan-Mayer ap- 
proach [ p^ , pT| in which effective potentials between pro- 
teins are invoked to explain crystallisation and the effects 
of added salt, often in the context of DLVO potentials 
JlOt where Debye-Hiickel theory is used to account for 
the charge interactions. However, a careful study of the 
relationship between the two approaches made by Hill in 
the 1950s indicates both approaches should converge on 
the same results in the high salt limit p^ . 

To apply these ideas to the protein crystallisation prob- 
lem, one constructs a free energy which enforces the 
charge neutrality constraint. Small ions can be treated as 
ideal solution species, but one must at least incorporate 
non-ideality of the protein solution, to allow it to form 
the ordered phase whch represents the protein crystals. 
The first thing to do is to set up this general theory. 



II. GENERAL THEORY 

I suppose the proteins under consideration have a 
molecular volume Vp, for example lysozyme is a charged 
globular protein of approximate size Vp = (tt/G) x 4.5 x 
3.0 X 3.0 nm'^ — 21.2 nm'^ [|lO|, or a molar volume 
NhVp = 12.8 M~^ where A^a is^vogadro's number. If 
the protein concentration is Cp, the effective volume frac- 
tion \s (f) = VpCp. In the absence of charge effects, a base- 
line model which incorporates a freezing transition to cor- 
respond to protein crystallisation will be fully specified 
by a dimensionless free energy density Z^"-* = V^/JF^^'/y, 
where F'^'^^ is the free energy of the protein solution or 
crystal, V is the system volume, and (3 = 1/kT. There 
may be different branches of /^"^ to correspond to the 
fluid and ordered (crystal) phases. Below I shall consider 
two possibilities for this baseline model: Sear's model of 
hard spheres with sticky patches, and the AHS model. 
The full free energy is obtained from the baseline model 
by adding in contributions from the colons and coun- 
terions (which I suppose to be univalent) and imposing 
charge neutrality. If the added salt concentration is Cg, 
there will be colons at a concentration Cs and counterions 
at a concentration Qcp + Cs- Denoting the dimensionless 
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free energy density for the full model by /, I write 
/(</), = + [log 

+ (g0 + 0,)[iog(Q0 + </),)-i] 

where the last two terms are the ideal mixing terms from 
the colons and counterions respectively, written using 
(j)s = VpCs for notational simplicity. Unimportant con- 
stants and terms linear in cj) or <j)g have been dropped. 
Non-ideality of the small ions is neglected in the present 
treatment, although this is certainly a refinement which 
should be considered for more quantitative work. 

This free energy is a function of two density variables: 
(j) and (j)s. A phase equilibrium such as between protein 
crystals and protein solution corresponds to equality of 
osmotic pressure and chemical potentials for both compo- 
nents. As a consequence there will in general be different 
values of (j)s in coexisting phases, a salt repartitioning 
effect which is not usually taken into account. To solve 
for phase equilibria, it is useful to transform / into a 
semi-grand potential ^i^"*) where ^i'^' = VpCs^'^ is a 
dimensionless salt reservoir concentration po|] . To make 
this transformation note that the salt chemical potential 
is 

df 

Pt^s^ j^^\og^s+\og{Q^ + (l)s). (2) 

ops 

Now (j)s (/>i^'' in the limit </) — > 0, therefore Pfi^ = 
2 log (^i^^ . This gives 

((/.f'f = (/.,(Q0 + (3) 

which is the usual Donnan equilibrium result that the 
product of the colon and counterion concentrations takes 
a constant value in all phases including the reservoir. 
Solving this gives 

0. = [(OV + 4(0f^)')'/'-Q0]/2. (4) 

The semi-grand potential h = f — I3^s4>s and the first 
two derivatives with respect to (p at constant (ps are, 
after a few lines of calculus, 

h = /(°^ + g01og(Q0 + Ps) - {Q<l> + Hs). (5) 

dh/d(i>^df^''yd^ + Q\og{Q(j) + c^s), (6) 

= 92/(°V50' + QVm + 20s). (7) 

The advantage of this transformation is that h is effec- 
tively a one-component free energy and can be treated 
accordingly. To use these results, one should remember 
to substitute for (ps from Eq. (|^). For example, the os- 
motic pressure follows from Vp(3Il = (t>{dh/d(j)) — ft,, ie 

n = n(°) +/cT(Qcp-f 2c,). (8) 

This shows explicitly the small ions behaving ideally, at 
a total concentration Qcp -\- ICs- 



What is obvious from these results is that there is a 
cross-over in behaviour at Qcp Cg or salt concentra- 
tions of the order Q(j)/Vp. The difference in protein vol- 
ume fraction between the solution and crystal is often 
A(/) - 0.5. Putting Q ~ 10, and Vp ^ 10 M^^, this corre- 
sponds to a cross-over salt concentration ^ 0.5 M. If the 
salt concentration is much less than this, there will be a 
large osmotic pressure difference between the crystal and 
the solution due to the counterions, having the effect of 
narrowing the coexistence gap. This is the basic reason 
why the coexistence boundary occurs at salt concentra- 
tions of this order of magnitude in this model. It is a 
much larger cross-over salt concentration than intuition 
might have suggested based on experiences for colloidal 
systems (eg Q - 10^, Vp ~ 10^ nm^ gives lO'^M), 
but then globular proteins are much smaller than colloids. 

On the other hand, typical protein solutions on the 
crystallisation boundary have (j) ~ 0.05, corresponding 
to a cross-over salt concentration ~ 50 mM. For salt 
concentrations much larger than this, the effects of salt 
and charge are subsumed into a scaling variable ^l^'/Q^. 

Since (jxs is more or less the salt concentration in the 
solution if the protein concentration is small, this is a 
natural explanation for the scaling properties described 
in the introduction. In the next sections, I place these 
arguments on a firm footing, starting with discussions of 
the second virial coefficient and the high-salt scaling of 
solubility, since these are not dependent on any particular 
baseline model. 



III. SECOND VIRIAL COEFFICIENT 

Any theory for the free energy of a protein solution 
contains a prediction for the second virial coefficient. 
For the present theory, B2 can be obtained from the 
osmotic pressure result Eq. (recalfing that Eq. (^) 
should be used for (/)„), or directly by inspection from 
Eq. (|). Either way, as advertised in the introduction, 
B2 = -\- /Acs where B^'^ is the contribution from 
the baseline model. One can prove the same result holds 
in light scattering determination since charge density 
fluctuations are suppressed in the long wavelength limit. 
We can hope that the baseline model is insensitive to 
the net charge and salt concentration, so that b'^^ is 
independent of Q and Cg, but is there any evidence for 
this? For lysozyme, Egelhaaf and Poon have collected 
experimental data on B2 from the available literature 
||l2|. This data is used to construct Fig. |l|(a) which 
shows B2 — Q^/Acs normalised by the value expected 
if lysozyme proteins behaved as hard spheres, namely 
B^s ^ _ 85 nm^ (this is close to Bl^^ = 82.3 nm^ 
used in Ref. ||2l|) Inspecting Fig. |l|(a), it appears there 
is indeed a data collapse to an approximate plateau for 
Cs 0.25 M, although a downwards trend of 40% or so 
can be detected at high salt and there may also be a 
problem if the protein charge is too small (the 'rogue' 
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FIG. 1: (a) sfV-Bl^^ for lysozyme, where = B2 

/4cg, as a function of salt concentration Cs [Q. The dotted 
line is the average, Eq. (^), for all data with Cs > 0.25 M. (b) 
Solubility data for lysozyme from Guo et al [11]. The dashed 
line is the least-squares fit to the data, Eq. (12) in the text. 



point with B^"^ > is for Q = 3.6). There is consid- 
erable variation in the data, even for measurements at 
supposedly identical state points, which reflects the ex- 
perimental difficulties in obtaining B2. Averaging over 
all results in the plateau region in Fig. ^a) arrives at 



B. 



(0) 



(-2.7 ±0.2) X B 



HS 
2 • 



(9) 



This result can be used as a constraint for any particular 
baseline model for lysozyme. Furthermore, the approxi- 
mate constancy of B^^ as Q and Cs are varied suggests 
that it is at least reasonable to make the baseline model 
independent of Q and Cg, as is assumed in the remainder 
of this paper. 



IV. HIGH SALT SCALING BEHAVIOUR 

In the high salt limit, a simple result can be obtained 
for the crystallisation boundary (equivalently the protein 
solubility) by assuming that the baseline model is inde- 
pendent of salt and charge, the density in the crystal 
remains constant, and the protein concentration in the 
solution phase becomes very small so that the solution 
can be treated as ideal. These features are certainly seen 
for the two specific baseline models discussed below. The 
simplifications amount to treating df^'^ydcj) in Eq. (||) as 
a constant independent of Q and Cg in the protein crys- 
tal, and taking the dilute limit df^^ydcj) — > log^p in the 
solution where (j)-p is the protein concentration in the so- 
lution phase. Equating protein chemical potentials in the 
two phases results in 



log (/)F + Q \og{Q(t>F + 0s, f) 

= constant 4- Q log((50x 



(10) 



where 0s^f and i/i^.x are the salt concentrations in the 
solution and crystal phases respectively, and 0x is the 
protein volume fraction in the crystal. The terms in this 
can be expanded assuming that (j)-p <^ 1, (ps^p « 0s^\ 
4's,x ^ Q(t>x, and making use of Eq. (H) to get 



log( 



constant 



(11) 



The prediction is that the logarithm of the protein solu- 
bility should be proportional to the square of the protein 
charge, and inversely proportional to the salt concentra- 
tion in the solution. Fig. |l](b) shows that such a law is 
indeed satisfied for the solubility data for lysozyme from 
Guo et al The best fit hue in Fig. 0(b) is 



/ solubility \ 

log( -r) 

^ mg ml ^ 



-0.25 ±0.15) 

+ (0.010 ±0.001) X 



(12) 



Cs/M 



The predicted slope from Eq. ( |ll| ) is (/)x/2T^. If we take 
0x ~ 0.5 as a reasonable estimate of the protein vol- 
ume fraction in the crystal, then (f)^/2Vp « 0.02 M which 
around twice the measured slope. However there are a 
number of omitted effects in the present theory which 
could account for the discrepancy. Nevertheless the ob- 
servation that the logarithm of the protein solubility is 
inversely proportional to the inverse salt concentration is 
apparently quite common . Also note that solubility 
~ indicates a very strong dependence on the protein 
charge. This is a natural explanation for the dramatic de- 
crease in solubility around the isoelectric point described 
in Ref. §. 



V. TWO SPECIFIC BASELINE MODELS 

For lower salt concentrations, the protein solubility is 
not small and the above treatment needs refinement. I 
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now consider the minimal extension to the basic model 
that takes into account protein non-ideality by specifying 
a particular baseline model. In this way, the advantage 
of an analytic transformation to the semi-grand poten- 
tial h in Eqs. is retained. The baseline model 
should encompass both the fluid and crystal phases, and 
have a second virial coefficient consistent with the anal- 
ysis above. One such suitable model has been devised by 
Sear. 



TABLE I: Fluid-fluid critical points for Sear's model showing 



critical parameter values (f)^'^'' , fc'*^-* and B. 
of Us- 



(C) 



for several values 



ris 


0(C) 


fc(C) 


Bf'/Bf' 


4 


0.09 


16.8 


-7.38 


5 


0.12 


7.02 


-3.39 


6 


0.15 


3.90 


-1.93 


8 


0.21 


1.77 


-0.77 



A. Sear's model 

Sear's model comprises hard spheres (HS) with 'sticky 
patches' The HS diameter is cr, chosen such that 

— Vp (ie cr ~ 3.4 nm for lysozyme). There are 
sticky patches per sphere which associate in pairs and 
are characterised by a range Tc > cr, an angular width 6c, 
and a depth e. The free energy of the fluid (F) phase (I 
reproduce only the essential details of the model here) is 

/^W = /hs(0) + 0n,[logp+ (1 -p)/2]. (13) 

In this /hs is the HS fluid free energy and p is the pro- 
portion of non-bonded sites, solving (1 — p)/p'^ = ^05hs 
where 



(14) 



is a dimensionless bond association constant. The asso- 
ciation equilibrium includes an enhancement factor, 5hSj 
for the HS pair correlation function at contact. The sec- 
ond virial coefficient in this model is given by 



r(°) - R 



HS 



- (n,/2) kVp. 



(15) 



It is apparent that the fluid phase properties are com- 
pletely determined by k and Ug. The model predicts 
fluid-fluid phase separation for sufficiently large values 
of fc, and Table | gives the critical point for ris = 4- 
8. Note that the second virial coefficient at the critical 
point in this model provides an marked counterexample 
to the observation of Vliegenthart and Lekkerkerker that 
B^2^^ / B^^ « —1.5 for a wide variety of other models [ p3[ . 
Whether fluid-fluid phase separation is metastable in this 
model depends on the actual values of rc, 9c and /?e, and 
examples of both are given by Sear. 

Sear provides a cell model for the protein chemical po- 
tential in an ordered cubic phase, arguing that osmotic 
pressure is not important. My approach here is slightly 
different. I use the same cell model to specify the free 
energy and take into account the osmotic pressure. The 
results are not essentially different from Sear's results, 
although the analysis does point up an important prop- 
erty of Sear's model (Fig. |^(a) below). Following Sear, 
the free energy per protein in the crystal (X) is 



/ a 



= C-31og - -1 

Vcr 



~2 



(3t w 



(16) 



where a/a — (6(/)/7r)~^/^ is the unit cell size relative to 
the HS diameter. The constant in this is C = log(lp/cr'^), 
provided that /hs ~ ^^(log 0—1) in the fluid phase as 
0^0. This free energy is appropriate for (j) ^ 4>min 
where ^min — 7r(cr/rc)'^/6 is the volume fraction around 
which the bonds in the crystal become dissociated. As 
(j) decreases past 0niin the last term in Eq. ( |l6|) vanishes 
rapidly. To implement this I have introduced an ad hoc 
cut-off function ii;[(0 — 4>u\in) / 5(1)] in the last term, where 
Scf) sets the rate at which the cut-off operates. For the 
present calculations I take w[x\ — 1/(1 + e~^) and Scf) = 
0.01. This cut-off function represents the way in which 
the short-range attraction potential falls off with distance 
in the model. The actual details may shift the phase 
boundaries but are not important for the broad picture. 

Although several values of the parameters in the model 
were examined, I only report in detail here on calcula- 
tions for Us = 6, Tc = 1.05cr and 9c = 0.45 (so the range of 
the attraction is rc—cr sa 2 A and the angular width about 
26°). These values were chosen to give quite good agree- 
ment with the crystallisation boundaries for lysozyme in 
the present model. Interestingly, in a separate analy- 
sis Curtis et al |2^ also conclude rig = 6-8 is appro- 
priate, and Oki et al identify three 'macrobonds' 
between lysozyme molecules from crystallographic data, 
again corresponding to ris = 6 contacts per protein. 

Fig. 11(a) shows the fluid and crystal free energies for 
Us = 6, = 1.05cr and 9c = 0.45, for two values of /3e. It 
is clear that there is a certain minimum value j3e ~ 4.8 
below which the crystal is metastable with respect the 
fluid. This is an important contrast with the AHS model 
(next subsection) in which the ordered phase is present 
at f3e — > 0. Crystallisation in Sear's model is essentially 
an energetic transition to an ordered phase dominated 
by the directional interactions (compare 'energetic fluid' 
concept introduced by Louis [p6|) and not a continuation 
of the entropic HS freezing transition. The position of 
the sticky patches controls the crystal structure, and as 
suggested by Sear, this may explain the relative ease or 
difficulty of crystallising various proteins. 

Fig. H(b) shows the phase behaviour for the model for 
the chosen parameter set. The solid lines in Fig. ^(b) 
are for the baseline model and include representative tie- 
lines. As pe increases, the fluid-crystal phase transition 
widens (the re-entrant fluid phase expected at larger 4> 
is not shown). There is a metastable fluid-fluid phase 




FIG. 2: (a) Free energies for baseline model illustrating 
metastability of crystal branch (X) with respect to fluid 
branch (F) for fie ^ 4.8. (b) Phase behaviour in (0, 13e) 
plane for full model: solid lines are for the baseline model 



{Q = ci^' — 0), dashed line is for Q 



10 and c 



(R) 



0, 



and chained lines are for Q = 10 and ci^' = 0.2 M. Other 



parameters are Ua = &, — 1.05cr, dc 



0.45. 



separation for fie > 7.18. 

The effect of charges and added salt is obtained by 
inserting the baseline model into the general formalism 
in section II. The dashed line in Fig. ||(b) shows the fluid- 
crystal phase boundary at Q = 10 in the absence of salt. 
The transition has been markedly narrowed and the fluid- 
fluid transition moves to such a high value of /3e that it is 
no longer on the diagram. Repeating the calculation for 
Cs^-* = 0.2M obtains the chained lines in Fig. §(b). The 
fluid-cystal phase boundary is intermediate between the 
zero salt limit and the baseline model, and the metastable 
fluid-fluid transition has reappeared in the diagram for 
/3e > 9.14. Finally as cl^' oo, the phase boundaries 
all move back to coincide with the baseline model. Thus 
the effect of charge in the model is to strongly suppress 
existing phase transitions in the absence of added salt. 



FIG. 3: (a) Phase behaviour in ((^, Cs) plane for Sear's model 
for Q = 10. (b) Fluid-crystal binodals for Q = 8.0 (triangles), 
9.4 (diamonds) and 11.4 (circles). Other parameters are Us = 
6, Tc = 1.05(7, dc = 0.45, f3e = 7.4. 



Adding salt weakens and eventually destroys the effect. 
The salt concentration required to do this is c., ~ Qcp, 
as discussed already in section II. 

I now set f3e = 7.4, marked by an arrow in Fig. ||(b), so 
—2.7 reproduces the value calculated in 
Fig. |^(a) shows a typical phase di- 



ih^iBf'>/Bf^ = 



[II above 



section 

agram for the full model at this value of /3e, as a function 
of added salt. Salt reservoir concentrations have been 
converted back to actual salt concentrations. The slop- 
ing tie-lines indicate a salt repartitioning effect in which 
the colon concentration is enhanced in the more dilute 
phases. In this simple model I have not taken into ac- 
count excluded volume in the dense phases, so in reality 
the salt repartitioning would be more marked (see below 
for a further discussion on this). In this representation 
we can clearly see how adding salt broadens the fluid- 
crystal coexistence. Metastable fluid-fluid phase separa- 
tion does not appear until salt concentrations ^ 2.3 M for 
this value of /3e. Note that the shape of the fluid-crystal 
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binodal (which I interpret to be the same thing as the 
crystaUisation boundary) is in agreement with common 
experience. This is a natural but non-trivial consequence 
of charge neutrality in the present model. 

This calculation is now repeated for Q = 8.0, 9.4 and 
11.4, which are the three values of lysozyme charge exam- 
ined experimentally by Poon et al. Fig. ^(b) shows the 
crystallisation boundary as a function of added salt, for 
the three values of Q. Clearly, the higher the charge, the 
more the phase transition is suppressed. Finally, these 
same phase boundaries are replotted in Fig. ^(a) as a 
function of the scaling variable Cg/Q^. In this represen- 
tation, the curves all collapse to lie on approximately the 
same quasi-universal crystallisation boundary. The scal- 
ing collapse is robust: if the calculations are repeated 
for different parameter values, the quasi-universal crys- 
tallisation curve moves up or down but a similar scaling 
collapse is always obtained. 

The mean crystallisation boundary from Fig. ^(a) is 
shown in Fig. |^(b), which now includes the experimental 
data from Poon et al. There is reasonable agreement on 
the location of the quasi-universal crystallisation bound- 
ary, although it appears there is always some discrep- 
ancy between the shape of the theoretical and experi- 
mental boundaries at low (j) which recalls the discrep- 
ancy in the slope of the high-salt law in Fig. ^b). It 
is worth re-emphasising that the present model is con- 
strained to have the correct i?2, at least for > 0.25 M. 
The absence of a metastable fluid-fluid phase separation 
agrees with Poon et al who observe that no such phase 
transition occurs at the temperature of the experiments 
(22.5 °C). However, a temperature decrease of just a few 
percent in the model, so that /3e = 7.7 for example, is 
sufficient to bring fluid-fluid phase separation to accessi- 
ble salt concentrations in the range 0.5-1.0 M. This is in 
accordance with the observations of Muschol and Roscn- 
berger |^ (but additional caution is required since it 
is unlikely the only effect of temperature is through the 
value of /3e). In the scaling representation of Fig. ^(a) 
such metastable fluid-fluid binodals also collapse to a 
quasi-universal curve — this is a prediction of the theory 
that would be interesting to test experimentally. 

The scaling collapse in Fig. ^ occurs because the ef- 
fects of charge and salt in fluid free energy are effectively 
combined into a single scaling variable Cs/Q^, at salt con- 
centrations in the vicinity of the crystallisation boundary. 
To show this is not just a feature of Sear's model, I now 
turn briefly to the more commonly studied AHS model. 



B. Adhesive hard sphere model 

The adhesive hard sphere (AHS) model is also a suit- 
able baseline model for the general theory in section II. 
In this model, hard spheres interact with a short range 
isotropic attractive potential. As noted by Rosenbaum et 
al 1^ , the phase behaviour is largely insensitive to the de- 
tails of the potential provided the second virial coefficient 




(b) t> 

FIG. 4: (a) Data of Fig. |(b) replotted in ((/>, Q"^Cs) plane, 
(b) Experimental data of Poon et al jlo|, for Q = 8.0 (trian- 
gles), 9.4 (diamonds) and 11.4 (circles). Filled (open) sym- 
bols indicate the occurence (non-occurence) of crystals. The 
dashed line is the approximate experimental crystallisation 
boundary. The solid line as the mean crystalisation bound- 
ary from (a). 

is used as the effective temperature axis. My approach to 
the baseline model here is closely based on that of Noro 
et al [|2^ who investigated the effects of long range forces 
on the AHS model. 

For the fluid phase, I use Barboy's treatment of Bax- 
ter's analytic theory Q. In Baxter's theory the 
attractive potential is characterised by a 'stickiness' pa- 
rameter T which is related to the second virial coefficient 
by 

i?f /i?2HS = l-(4r)-i. (17) 

The crystal phase is expected to be FCC since this has 
the greatest density of intersphere contacts. As the stick- 
iness is switched off (r ^ oo or (3e 0) the fluid-crystal 
phase transition goes over into the usual HS freezing tran- 
sition, which is a contrast to Sear's model. The AHS 
ordered phase has always proved rather more difficult to 
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FIG. 5: (a) Phase behaviour in {(p, r) plane for the adhesive 
hard sphere model with n = 30: solid lines are for baseline 
model with Q = ci^' — 0, dashed line is for Q = 10 and ci^' — 
0. The shaded region is the 'danger zone' Eq. (|l9|). (b) Fluid- 
crystal binodals in {(l),Q~^Cs) plane for Q = 8.0 (triangles), 
9.4 (diamonds) and 11.4 (circles); other parameters are n — 
30, r = 0.1. 



treat analytically than the fluid phase, and approaches 
have ranged through density functional calculations by 
Marr and Gast Q, and Tejero and Bans Q, to detailed 
simulation studies . Here I use the cell model of Daa- 
noun et al p2| for the crystal phase free energy, assuming 
the attractive part of the AHS potential is — e(r/cr)^" 
with n 3> 1. The cell model free energy is almost iden- 
tical to that already used in Eq. (|l^ ) with = 12 for 
the FCC structure. The difference is that there is no 
restriction on orientation so the term — log(0j?/7r^) is ab- 
sent, and the role played by the ad hoc cut-off function 
is taken over by the functional form of the actual short 
range potential, ie w[x] {a/a)^". Like the choices for 
the cut-off function in the previous section, the actual 
value of n may shift the phase boundaries but does not 
change the broad picture. For the purposes of the present 
calculation I set n = 30. 



A link between the cell model for the crystal phase and 
the Baxter-Barboy treatment of the fluid phase can be 
made by matching B2 of Eq. (|^ with the exact expres- 
sion 

b(o)/BHS ^ 1 + 3 / dxx^[l- exp(-/3ea;~")]. (18) 



Fig. ||(a) shows the phase behaviour for the model. The 
solid line in this figure is the fluid-crystal phase boundary 
for the baseline model, and the dashed line is the corre- 
sponding boundary for Q — 10 with no salt. Similar 
to Fig. H(b), the effect of charge is to narrow the fluid- 
crystal transition. As salt is added (not shown here), the 
boundary moves back towards that of the baseline model. 

There is a complication that arises for Baxter's solu- 
tion for the AHS fluid phase free energy. Baxter's theory 
involves the solution of a quadratic equation. If 



T < 



(12(/. + 60^)1/' -60 
6(l-(/.) 



(19) 



then this quadratic equation has complex roots and Bax- 
ter's theory becomes inadmissable. Normally this 'dan- 
ger zone', shown shaded in Fig. ^(a), is happily hidden 
within the fluid-crystal two phase region. In the charged 
version though, it emerges into the single phase fluid re- 
gion. For this reason one cannot choose r < (2 — •^2)76 « 
0.098, the maximum of Eq. (l9h, and unfortunately this 



excludes r = 0.067 which would match B^2^ /B. 



from section HI 



'2 



-2.7 



To demonstrate the scaling collapse for this model, I 
therefore choose r = 0.1 (equivalent to B^"^ /B^^ — —1.5 
psf.) I repeat the calculations of the previous section to 
obtain the crystallisation boundaries (fluid-crystal bin- 
odal) as a function of added salt, for Q = 8.0, 9.4 and 
11.4 as used previously. Fig. ||(b) shows these bound- 
aries plotted using the scaling variable c^/Q^. Again the 
curves collapse to a quasi-universal crystallisation bound- 
ary, similar to those in Fig. |^(a). This result confirms 
that the scaling collapse is a common feature for two dif- 
ferent baseline models for the phase behaviour. 



VI. SALT REPARTITIONING AND 
INTERFACE STRUCTURE 

These ideas have some additional interesting physical 
implications for the equilibrium between crystals and so- 
lution which are explored in this section. The first impli- 
cation concerns salt repartitioning. In the simple theory 
of section II, the product of colon and counterion concen- 
trations is a constant as in Eq. (||) . Each phase has to be 
electrically neutral, so that a phase enriched in protein is 
also enriched in counterions, and consequently depleted 
in colons (Donnan common ion effect). This explains the 
slope of the tie-lines in Fig. ||(a) for example. We can ask: 
is there any independent evidence for this phenomenon? 
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TABLE 11: Lysozyme salt repartitioning data from Palmer 
et al H. Solution is 0.855 M NaCl / 0.2 M NaAc buffer 
{pH = 4.5). Ion concentrations in the third column are calcu- 
lated by multiplying those in the second column by 1/(1 — 0), 
where (j) = 0.64 is the crystal volume fraction. The charge 
unaccounted for in the crystal, assuming Q = 11.4, is 0.05 M 
or < 3% of the total charge. 



cone / M 


solution 


crystal 


crystal free volume 


Lysozyme'^"'" 




0.056 




Ac" 


0.2 






Na+ 


1.055 


0.19 


0.53 


Cl~ 


0.855 


0.78 


2.2 


Na+ X Cr 


0.9 




1.2 



Let me preface the answer to this by some cautionary re- 
marks. Deviations from the simple Donnan equilibrium 
result may arise from three sources. Firstly excluded vol- 
ume effects mean one should really consider the small ion 
concentrations in the available free volume. This will be 
an important consideration in the crystal phase where 
the volume fraction occupied by the protein is ~ 50%. 
Secondly, deviations from ideal solution behaviour can 
be expected at high salt concentrations — such deviations 
are usually absorbed into 'activity coefficients'. Thirdly 
there may be significant specific ion effects, such as seen 
in the Hoffmeister series. 

Repartitioning of salt was studied in detail by Vekilov 
et al [Q. They found non- uniform salt concentrations 
inside protein crystals, but consider this is likely to be an 
effect of impurities or growth kinetics. The present the- 
ory only addresses the equilibrium repartitioning of salt 
in perfect lysozyme crystals though. Fortunately, some 
relevant equilibrium results for salt repartitioning can be 
found in recent experiments by Morozova et al [ p5[ on 
cross-linked lysozyme crystals in contact with salt solu- 
tions. For example. Table || shows concentrations of Na+ 
and Cl~ in solution and in the crystal from the early work 
of Palmer et al cited by Morozova et al. There is 
marked repartitioning, for example the colon concentra- 
tion in the crystal free volume is about half that in the 
external solution. The product of the two ion concen- 
trations is approximately constant though (final row in 
Table ^ , provided the protein excluded volume is taken 
into account. Detailed calculations by Morozova et al 
for their more recent data also take into account activity 
coefficients and give excellent agreement between theory 
and experimental results: they conclude ". . . for small 
ions capable of penetrating into the crystal channels [the] 
electrostatic (Donnan) potential controls the equilibrium 
internal concentration of ions in just the same way as in 
polyelectrolyte gels" [Q. By way of contrast, the same 
work also suggests Br~ has a significant specific interac- 
tion with lysozyme. 

The second implication concerns the electrical struc- 
ture of the interface between crystal and solution phases. 



In a Donnan membrane equilibrium [ p^ [TJ, , an elec- 
trostatic potential difference, the Donnan potential, de- 
velops between the two compartments on either side of 
the membrane. In the present theory for protein crys- 
tallisation, it is still true that a difference in the mean 
electrostatic potential develops between the crystals and 
the solution phase. Whether this should also be called a 
Donnan potential, a Galvani potential, or perhaps some- 
thing else, can be debated p7[ |. 

At any rate, the potential difference is readily calcu- 
lated given the salt repartitioning. If c+ and c_ are the 
colon and counterion concentrations respectively, then 
c± = ci^^ exp[^/3e^], where Jp is the mean electrostatic 
potential in the phase of interest measured relative to 
the salt reservoir (this result can also be used to re- 
cover Eq. (^). The difference in mean potential between 
fluid (F) and crystal (X) phases, Aip = Tp^^^ —lp^^\ is 
A(/5 — {j3e)^^ \og{cf^ I c's"''). For example, for the salt 
concentrations in Table |l| cf ^ = 1.055 M, ci^^ = 0.53 M 
(in the free volume), and therefore A(/9 « 17 mV. The 
potential difference arises because an electrical double 
layer is formed at the crystal-solution interface and is in- 
timately connected with salt repartitioning. The details 
of this will be discussed in a forthcoming paper |Q . An- 
other consequence of the Donnan potential is that the pH 
in the crystal will be 0.434A(y9 higher than in the solu- 
tion (where 0.434 = logj^g e). If the solution pH is below 
the isoelectric point, as is often the case for lysozyme, 
one might expect this to lower the charge per protein in 
the crystal although it is unlikely that the charge can be 
determined with any accuracy from such a naive calcu- 
lation. 

Pair-potential theories miss both the effects of salt 
repartitioning and the potential difference between the 
two phases, which are essentially many-body phenom- 
ena H, 0). Whilst this is not a problem at high 
salt where Hill's mapping between the McMillan-Mayer 
/ Debye-Hiickel pair potential approach and Donnan's 
method goes through |l^, the effects can be significant 
at lower salt concentrations. The importance of the ef- 
fects can be gauged by comparing the magnitude of At/? 
to the thermal energy kT. If A(/9 > kT /e then there 
may be significant errors introduced by using a single 
effective pair potential (such as a DLVO potential) for 
both the crystal and the solution phases. On the other 
hand, if A(^ ^ /cT/e, the use of pair potentials cannot 
be criticised on these grounds. It can easily be shown 
that A(/3 ^ kT /e corresponds to the by now familiar 
Cs ^ Qcp, in other words salt concentrations smaller than 
the protein charge density in the crystal, which is typi- 
cally 0.5 M. 



VII. CONCLUSIONS 

There are several conclusions from the present work. 
The first concerns the lysozyme / NaCl system. The 
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quasi-universal crystallisation boundary observed by 
Poon et al can be fit reasonable well by the present 
version of Sear's model with ^ 6 sticky contacts per 
molecule, similar to the conclusions of other workers 
p^ , [2^ . The scaling collapse for B2 noticed by Egel- 
haaf and Poon is here attributed to a contribution 

Q^/4cs added to a 'bare' Bj"'' which is approximately in- 
dependent of salt concentration above about 0.25 M. In 
this system therefore, NaCl appears to be acting as an in- 
different electrolyte in the sense that it does not seem to 
exhibit specific ion effects. This appears to be confirmed 
by the experiments of Morozova et al discussed in 
the preceeding section. Other salts, for example NaBr, 
may not behave the same way of course. 

The other conclusions are general ones. Firstly, a sim- 
ple extension of existing models to incorporate salt ions 
and charge neutrality provides a straightforward expla- 
nation for the shape of protein crystallisation boundaries 
and the associated scaling properties seen for lysozymc. 
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